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We study the A-nucleon (AA^) force by using lattice QCD. The Bethe-Salpeter ampUtude is cal- 
culated for the lowest scattering state of the AA^ so as to obtain the AN potential. The numerical 
calculation is twofold: (i) FuU lattice QCD by using 2 + 1 flavor PACS-CS gauge configurations 
with, /3 = 1.9, corresponding to the lattice spacing of a = 0.0907(13) fm, on a 32^ x 64 lattice. 
A set of parameter (K;,f/, JQ.) = (0.13770,0.13640) is used, which corresponds to m„ « 300 MeV 
and niK ~ 594 MeV. The spatial lattice volume corresponds to (2.86 fm)^. (ii) Quenched lattice 
QCD with j3 = 5.7, the lattice spacing of a = 0.1416(9) fm, on the 32^ x 48 lattice. Two sets 
of hopping parameters {k„ii, k^) — (0.1665,0.1643), (0.1670,0.1643) are used. The spatial lattice 
volume is (4.5 fm)^. For the full QCD, we find that the Ap has a relatively strong (weak) repulsive 
core in the ^So (^Si) channel at short distance, while the potential has slight attractive region at 
medium distance. The lowest scattering energy in the finite lattice volume is calculated; Slightly 
negative values obtained in both spin channels. For the quenched QCD, we find that the results 
are qualitatively in agreement with those in the full QCD calculation. 
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1. Introduction 

Study of the hyperon-nucleon (YN) and hyperon-hyperon (YY) interactions is one of the im- 
portant tasks in the contemporary nuclear physics. These interactions are the bases to explore the 
strange nuclear systems, in which hyperons (or strange quarks) are embedded in normal nuclear 
systems as "impurities" Various issues, such as spin- and flavor-dependence, effect of the fla- 
vor SU (3) breaking, effect of hyperonic mixing (e.g., AN — LN and AA — EN), are the topics to 
be solved experimentally and/or theoretically. For example, recent systematic study (e.g., Ref. [^) 
for light (i'-shell) A hypernuclei (^H, ^H, ^He and ^He) suggests that the AN interaction in the 

channel is more attractive than that in the channel. However, the present YN and the YY 
interactions have large uncertainties, because the scattering experiments are either difficult or im- 
possible due to the short life-time of hyperons. Theoretically, it should be desirable to understand 
the YN and YY interaction (or, in more general, baryon-baryon interaction) based on the dynamics 
of quarks and gluons as fundamental degrees of freedom. If one can perform such an appropri- 
ate deduction along the theory of quantum chromodynamics (QCD), they should have a reliable 
prediction regarding the YN and YY potentials. 

Recently, the lattice QCD studies have been performed for not only the NN potential but also 
the EN potential in quenched QCD [^, ^]. See also Refs. ^ ^, 0] for the recent developments. The 
lattice QCD calculations focusing on the scattering parameters for AN system based on Liischer's 
formula have been reported in Refs. ||8|, ^]. The purpose of this report is to calculate the AN 
potentials by using the full and quenched QCD gauge configurations. The main results of this work 
are obtained by using the 2+1 flavor PACS-CS gauge configurations with the spatial lattice volume 
(2.86fm)^ jro[]. We also report the results calculated by using quenched lattice QCD with larger 
spatial lattice volume (4.5fm)^. 



2. Formulation 



The basic formulation is already given in Refs. [|3[ |[|]. (See also Refs. [0 [l|].) In this 
report, we briefly explain the procedure to obtain the potential from lattice QCD together with 
some technical point newly introduced in the present work. We start from an effective Schrodinger 
equation for the equal-time Bethe-Salpeter (BS) wave function: 



(2.1) 



Here ji. = m/^mN / {m/^ + miq) and E = /{2}Jl) are the reduced mass of the AN system and the 
non-relativistic energy in the center-of-mass frame, respectively. We consider the low-energy scat- 



tering state so that the nonlocal potential can be rewritten by derivative expansion [13], U(f,r' 



VANif^:^)5{r — r'). The general expression of the potential Van is known to be [14] 



Van = Voir) + V„ir){dA-dN) + VT{r)Sn + VLsir){L-S+) + VALs{r)iL-S-) + 0{V^). (2.2) 



Here Su = 3(aA • n){dN -n) — Oa-On is the tensor operator with n = r/\r\, S± = {On ± aA)/2 ai^e 
symmetric (+) and antisymmetric (— ) spin operators, L = —if x V is the orbital angular momentum 
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operator. We note that the antisymmetric spin-orbit forces (Vals and Valsx) do not arise in the NN 
case because of the identical nature of the nucleon within the isospin symmetry. 
First we calculate the four-point correlator by using the lattice QCD 



0), (2.3) 



where the summation overX is to select the state with zero total momentum. The Pa{x) and Ap{y) 
denote the interpolating fields for proton and A 

Paix) = eabc{Ua{x)C'Y5db{x))Uca{x), (2.4) 

^fi{y) = £abc {{daiy)Cr5Sh(y))uci}{y) + {sa{y)CY5Uh{y))d^p{y) -2{ua{y)CY5db{y))s,.p 
and A' pi {to) and p'a'{to) denote the wall source 

P'a'it) = X ^abc{Ua{xi,t)CY5dtix2,t))Uca'(x3,t), (2.6) 

^'p'it) = L ^abc {{da{h,t)CY5Sb{h,t))ucj}'{y3,t) + {sa{yut)CY5Ub{h,t))dcp>{y3,t) 

-2{u„{yy,t)CY5db{y2,t)Kp'{y3,t)} ■ (2.7) 



yi,y2,y3 



The Fast Fourier Transformation (FFT) algorithm is imposed into the actual C++ code to reduce 
the computational time, which makes the calculation at all of the spatial points (32^) possible. (In 
the previous work for the E^p system, we calculate only on the x-, y- and z-axes and the nearest 
neighbors in the long distance region. [§]) 

In the Monte Carlo calculations, noise reductions are made for the four-point correlator (2.3) 
obtained from lattice QCD in order to restore (i) the rotational (cubic group) symmetry, (ii) the 
spatial reflection symmetry, (iii) the charge conjugation and time-reversal symmetry. Then, we 
perform the spin projection for the source to obtain the wave function with particular spin compo- 
nent (7 =0,1). 

Fp,Ap{7,t -to;J,M) = Pa,p}Gap;a'P'{7,t -to), (2.8) 

a'j3' 

See Ref. [Q] for detail. Finally, the desirable wave function ^^.p is obtained from the projected 
four-point correlator at large t — to'. 



Fp^A^ir,t-to;J,M) = £a„£(0 pcc{X + r,t)ApiX,t 



X 



E„)e-^"^'-">'> (2.9) 



« Ao0„/3(r;/,M)e-^<'('-'») {t-to^ large). (2.10) 

Here £„ (\En)) is the eigen-energy (eigen-state) of the six-quark system with the particular quantum 
number {J^,M), and A„ = Y,a' P' Paip ^Fn\A! pi p' ■ For simplicity, we consider the effective 
central potential, Vc{r\J = 0) = Vo(r) - 3Va(r), for ^So and, Vc{r,J = 1) = Vo{r) +Vcy{r) + ■■■, 
for ^S\, where • • • is the higher order contribution from noncentral force such as Vjir), which is 
expected to be small in the present calculation. The effective central potential Vc{r;J) is obtained 
by focusing only on the S'-wave component of the wave function in each spin channel /: 

, , 1 V2(/)(r;7) 
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3. Numerical calculation 



3.1 A^/ = 2+l QCD 

Main results in this report are obtained by using the gauge configuration generated by PACS- 
CS collaboration with 2+1 flavor full QCD [jl^ ; The Iwasaki gauge action and the nonperturba- 
tively C?(fl') -improved Wilson quark action are employed. Calculations are carried out at j8 = 1.9 
on a 32^ x 64 lattice, corresponding to a = 0.0907(13) fm where a is the lattice spacing at the 
physical point | IC ] . A set of hopping parameters ( Kud , K^s 

) = (0.13770,0.13640) is used in this 
work. Several light hadron masses obtained by PACS-CS are shown in Table |I} The BS wave func- 
tion is obtained by employing the wall source at the time-slice to = 0, and the Dirichlet boundary 
condition is imposed along the temporal direction on the time-slice t = 32 with the Coulomb gauge 
fixing. The present result of the BS wave function is calculated with A'conf = 422, where A'^conf is 
the number of gauge configurations. 



3.2 Quenched QCD with larger spatial volume 

In this case, we use the plaquette gauge action and the Wilson fermion action with the gauge 
coupling j8 = 5.7 on the 32^ x 48 lattice. The periodic boundary condition is imposed for quarks 
in the spatial direction. The wall source is placed at = with the Coulomb gauge fixing and the 
Dirichlet boundary condition is imposed at f = 24 in the temporal direction. These setup are es- 
sentially the same as our previous calculations for the A^A^ and EN potentials Q] except for 
the treatment of the temporal part. The lattice spacing at the physical point is determined as 
a = 0.1416(9) fm {l/a = 1.393(9) GeV) from nip = 770 MeV. The hopping parameter for the 
strange quark mass is given by = 0.16432(6) from mx = 494 MeV. The spatial lattice volume 
is (4.5fm)^, which is enough to accommodate two baryons. The present result of the BS wave 
function in quenched QCD is obtained with A'^conf = 550. 



4. Results and Discussion 



4.1 Nf = 2 + \ QCD 

The left panel of Figure |l| shows the wave function obtained at the time slice t — to = 7. The 
red circle (blue triangle) corresponds to the ^So CSi) channel. In order to find the ground state 



Ma 



mi 



2+1 flavor QCD by PACS-CS with {K^d, 

Ref.rfM] 296(3) 848(20) 594(2) 



= (0.13770,0.13640) 

985(8) 1093(19) 1254(14) 



1315(15) 1448(10) 



quenched QCD with j8 = 5.7, k, = 0.1643, 

Kud = 0.1665 514(1) 859(3) 607(1) 904(3) 1287(8) 1345(6) 1368(6) 1412(5) 
K-„^ = 0.1670 465(1) 840(6) 588(1) 894(3) 1234(11) 1308(7) 1341(8) 1396(5) 



Exp. 



135 



770 



494 



892 



940 



1116 



1190 



Table 1: Hadron masses in the unit of MeV. The results for 2 
Ref. 0. A^conf = 560 in quenched QCD. 



1320 



1 flavor QCD by PACS-CS is taken from 
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Figure 1: Left: The radial wave function of Ap, in ^So (red circle) and "^^i (blue triangle) channels, obtained 
by using the PACSCS gauge configurations at t — to — 1 ■ A^conf ~ 422 gauge configurations are used. Right: 
The effective central potential for Ap, in the '^o (red circle) and (blue triangle), obtained from the wave 
function at time slice t — to — 7- The left (right) inset shows its enlargement in the ^Sq (^Si) channel. 



saturation of the Ap system, we define the "effective mass" of the wave function: 

m,n{t-to,T) ^ log ( , ) ■ (4.1) 
\FpA[r, {t + l)-to)J 

Figure ^ shows the effective mass at several spatial points. The plateaux starting appears at the time 
slice t — to = 6 in the channel (in the right panel in Fig. ^) while the plateaux starting seems 
to appear at the time slice around t — to = 1 in the ^5o channel (in the left panel). We need more 
statistics. In this report, we show the results obtained at the time slice t — to = 1. 

As is shown in Ref. [12], the non-relativistic energy E = k^/ (2/x) is determined by fitting the 
wave function in the asymptotic region in terms of the Green's function 

Gir^) = ^ I ^2^e'>"^ r = |p; p = G Z^l , (4.2) 

which is the solution of [A + k^)G{r,k^) = —8i{f) with 5f,(r) being the periodic delta function 
[11, 12]. We attempted to make a fit in the asymptotic region (15 < \r\ < 16) at the time slice 
t — to = l. The fitting region is chosen so that the potential becomes zero within the errorbars. The 
energy obtained in the ^So (^50 channel is E = -1.9(1.3) MeV (£ = -1.3(1.2) MeV). The present 
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Figure 2: The effective mass of the wave function for Ap at several spatial points {x,y,z), in the ^So (left) 
and the (right) channels, obtained by using the PACS-CS gauge configurations. A^conf = 422 is used. 
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results for the energies suggest that the spin dependence of the interaction in the low-energy region 
seems to be weak, although we need more statistics (and perhaps larger spatial volume) to make 
definite conclusion. The right panel of Figure |l| shows the effective central potential of Ap system 
obtained at the time slice t — tQ = 7. The strong repulsive core is obtained in the channel, while 
the repulsive core in the ^Si channel is relatively weak. The left (right) inset shows the enlargement 
in the (^S\} channel; A slight attractive well can be seen in the both spin channels. 

4.2 Quenched QCD 

We calculate the wave functions for two sets of the hopping parameters {Kud, K^) = (0.1665,0.1643) 
and (0.1670,0.1643). Table[I]also lists the hadron masses obtained from these quenched QCD cal- 
culations. The Figure || shows the wave function obtained with {Kud,Ks) = (0.1670,0.1643) at 
the time slice t — to = 7, by using Nconf = 550 gauge configurations. For both cases of calcula- 
tions with (K:„d, K,) = (0.1670,0.1643) and (0.1665,0.1643) in the and in the channels, the 
ground state saturations are achieved at the time slice t — to = 7. The energies obtained by fitting 
the wave function in the asymptotic region (10 < |r| < 16) to the Green's function (^.2) in the ^Sq 



CSi) channel are £ = -0.37(12)MeV (E = -0.47(ll)MeV) for (k-,,^, K-,)=(0.1670,0.1643), and 
E = -0.31(9)MeV {E = -0.38(7)MeV) for (k-„j, k;,)=(0. 1665,0. 1643), respectively. The lowest 
scattering energies in both and channels become lower as the u,d quark mass decreases. 
We hardly find the clear conclusion whether the interaction in the ^Sq channel is more attractive 
than that in the channel. Figure ^ shows the effective central potentials for Ap in the ^So CSi) 
channel in the left (right) panel, obtained from the wave function. For both (^^o and ^^i) spin 
channels, the height of the repulsive core increases as the u,d quark mass decreases. 

In summary, we calculate the AA'^ potential from lattice QCD. The full lattice QCD calculation 
for the AN system with larger volume and smaller lattice spacing at physical quark mass is highly 
desirable to clarify the spin dependence of the interaction. 
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Figure 3: Left: The radial wave function of Ap, in ^^o (red circle) and (blue triangle) channels, obtained 
by quenched QCD calculation with (jc,„/, K,) = (0. 1670, 0. 1643), at the time slice f — fo = 7. 
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Figure 4: Left: The effective central potential for Ap in the channel at m-,i ~ 465 MeV (red circle) and 
ffiji ~ 514 MeV (green box). Right: The effective central potential for Ap in the channel at m-,i ~ 465 
MeV (blue triangle) and m^i ~ 514 MeV (magenta diamond). The inset shows its enlargement. 
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